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We analyze the motion of a single hole on a Neel background, neglecting spin 
fluctuations. Brinkman and Rice studied this problem on a cubic lattice, introducing 
the retraceable-path approximation for the hole Green's function, exact in a one- 
dimensional lattice. Metzner et al. showed that the approximation also becomes 
exact in the infinite-dimensional limit. We introduce a new approach to this problem 
by resumming the Nagaoka expansion of the propagator in terms of non-retraceable 
skeleton-paths dressed by retraceable-path insertions. This resummation opens the 
way to an almost quantitative solution of the problem in all dimensions and, in par- 
ticular sheds new light on the question of the position of the band-edges. We studied 
the motion of the hole on a double chain and a square lattice, for which devia- 
tions from the retraceable-path approximation are expected to be most pronounced. 
The density of states is mostly adequately accounted for by the retraceable-path 
approximation. Our band-edge determination points towards an absence of band 
tails extending to the Nagaoka energy in the spectrums of the double chain and the 
square lattice. We also evaluated the spectral density and the self-energy, exhibiting 
k-dependence due to finite dimensionality. We find good agreement with recent nu- 
merical results obtained by Sorella et al. with the Lanczos spectra decoding method. 
The method we employ enables us to identify the hole paths which are responsible 
for the various features present in the density of states and the spectral density. 
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I. INTRODUCTION 



The discovery of high-temperature superconductivity has intensified the interest in the 
study of strongly correlated electron systems. Among the important related subjects await- 
ing solution is the dynamics of holes in the presence of a spin background. Among pioneering 
works on the problem one can cite that by Bulaevskii, Nagaev and Khomskii who, study- 
ing the motion of a hole on an antiferromagnetic background, with the neglect of quantum 
fluctuations, introduced the "string potential" concept. As the hole moves on the Neel 
background it leaves behind a string of overturned spins, thus increasing the exchange en- 
ergy proportionally to the length of the string. This can lead to self-trapped (or, as they 
were called in Ref. 0, "quasi-oscillator" ) states centered at the original position of the 
hole. Two years later, Brinkman and Rice studied the motion of a single hole on differ- 
ent spin backgrounds, employing the Nagaoka expansion for the hole propagator ||. This 
involves contributions by hole paths restoring the spin background, classified according to 
their length. For a cubic lattice and ferro-, antiferromagnetic and random spin backgrounds, 
they calculated the exact contributions of hole paths up to length 10, neglecting spin fluctu- 
ations. The density of states obtained was compared with the one resulting from considering 
only the self-retracing paths of the hole, the "retraceable-path approximation" (denoted RPA 
throughout this paper). In this approximation, which does not distinguish between different 
spin backgrounds, one can exactly sum the series for the propagator, which involves only 
local contributions, thus obtaining no momentum(k)-dependence: the motion is completely 
incoherent. With the RPA Brinkman and Rice H were able to account for the body of the 
band in the case of a Neel background. Less agreement was obtained || for the random 
background, while the ferromagnetic case reduces to the uncorrelated single particle motion 
in a tight-binding band which generally is quite different from the RPA description. They 
also observed that the RPA became exact in the one-dimensional case where, for any spin 
background, the free-particle result is obtained. The band edges appear at the so called 
Nagaoka energies. 
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In a series of recent works H|| the motion of a hole on a spin background was analyzed 
in the case of infinite dimensions. In particular, Metzner et al. ||] estimated the contri- 
butions to the propagator in terms of the dimension(D). Neglecting spin fluctuations, the 
only non-retraceable contributions which would remain finite at D=oo would be those of 
loop paths circulated only once. These would certainly not preserve a Neel ordered spin 
background, so that the RPA becomes the exact solution of the problem in this case. For 
the Neel background at finite dimensions, the lowest order corrections to the RPA density 
of states in the case of a hypercubic lattice would come from circulating elementary square 
plaquettes three times (thus of order ~ l/D 4 ). More important are the non local correc- 
tions to the propagator, the lowest of which is of order ~ l/D (the hole can propagate to 
nearest neighbour sites on the same sublattice circulating an elementary square plaquette 
one and a half times). The latter corrections would give rise to a k-dependent self-energy 
for the hole. Trugman || already mentioned these non-retraceable paths in his study of 
the two-dimensional Hubbard antiferro magnet. Quite different physics is expected when 
the quantum fluctuations associated with transverse exchange interactions are taken into 
account from the start. This problem has been addressed in many recent works, among 
others see Refs. [pH^C- 

Taking into account the results mentioned above for hole motion on a Neel background 
in the absence of quantum fluctuations, we propose here a new method to evaluate the 
local and non-local parts of the propagator in order to analyze the departures from the 
RPA for reduced dimensionality. In Section |I| we describe the method employed to obtain 
these quantities. Basically, the Nagaoka expansion of the propagator is rewritten in terms 
of non-retraceable skeleton-paths dressed by retraceable paths. Though one can explicitly 
evaluate the contribution of a dressed skeleton diagram of given length, one has to determine 
numerically the numbers of bare background-restoring skeleton-paths. The results obtained 
with this method for the double chain and the square lattice are presented in Section [[II]. 
We can quite reliably determine the band edges, and our results place them at an inter- 
mediate value between the RPA edge and the Nagaoka energy. In general, we find good 
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agreement with the recent numerical results of Sorella et al. fll4| , obtained with the Lanczos 
spectra decoding method. The method we employ enables us to determine the origin of the 
distinctive features exhibited by the density of states and the spectral density. Thus, we can 
identify the relevant paths responsible for the k-dependence at reduced dimensionality. In 
Section [TV] we summarize our results. 

II. ANALYTICAL METHOD 

In this section we describe the method employed to evaluate the hole propagator on a 
Neel background, in which spin fluctuations are neglected. 

The Nagaoka expansion |3|,|4j] for the Green's function G describing the propagation of a 
hole inserted into an arbitrary fixed spin configuration | s >, takes the following form: 

G U W ) = E < s I c t-^-FT^ I s > 



1 

w 



l ' a w-H 

n=l v w 



H = -t Yl [(! -^-^taCjA 1 - n 3,-*) 
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+ h.c.]. (1) 

Above, cf a (q j(J ) creates (annihilates) a particle with spin o at the lattice site i, while H 
is the Hamiltonian describing the correlated hopping of a particle between nearest neighbour 
sites. The coefficients Af ■ n denote the number of distinct n-step paths of the hole in the spin 
background s which start at site j and end at site i, restoring the original spin configuration 
after the last step. In particular we will consider the case of a Neel spin background, s = N, 
where propagation is only possible between sites belonging to the same sublattice. 

As we mentioned in the previous section, the paths in which the hole completely retraces 
its steps are the only relevant ones at D = 00 on a Neel background M. Even at finite 
dimensions they would determine the leading terms of the local Green's function. Following 
Brinkman and Rice P| one can derive the RPA propagator by introducing the following 
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irreducible "self-energy" S (here irreducibility means that the path cannot be split up into 
two or more consecutive retraceable paths): 

1 



G RPA (w) 



(2) 



w(l-S(w)) 

The lowest order RPA self-energy will derive from the contributions of single jumps to a 
nearest neighbour and immediate return to the origin, so that: 



S°(w) 



w 



(3) 



where Z is the number of nearest neighbours. The next order irreducible self-energy con- 
tribution consists of one jump to a nearest neighbour of the origin followed by a further 
jump in any direction different from the previous jump (thus having [Z — 1) possibilities for 
this step) before retracing the whole path. In general, one can write the RPA irreducible 
self-energy as: 

t 



s(w) = z(-yc(w), 

w 

C(w) = i + (^Z±)S(w) 
1 



1 



z 



) 2 S(w) 2 + ... 



(4) 



(5) 



l-[^)S{w) 

that is, it is written self-consistently in terms of the jump to a nearest neighbour of the 
origin dressed by retraceable paths, where C is the RPA-dressed irreducible vertex part. 



The solution of the self-consistent system of equations (|4]) and 

Z 



is 



S(w) 



2(Z - 1 

Therefore the RPA propagator is given by: 
G RPA {w) = 



1-Jl-4(Z-1)(^ 



(6) 
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From this expression one can easily determine the band edges, ±wrpa, of the energy 
spectrum obtained in the RPA for a hole hopping on a Neel background as: 
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wrpa = 2VZ - It, 



(8) 



which exhibits a reduction from the free-particle value (Nagaoka energy): wpp = Zt for 
Z > 2, while the RPA coincides with the free-particle result in the one-dimensional (Z = 2) 
case, as already mentioned in the previous section. 

Having introduced the retraceable-path approximation (RPA) which gives a good ap- 
proach to the local propagator, we will now describe the method we employed to obtain the 
Green's function of a hole moving on a Neel background. In the Nagaoka expansion ([[]) we 
take into account all background-restoring paths in the following way: for each path-length, 
we take all possible non-retraceable "skeleton" -paths which restore the Neel background, and 
dress them by inserting all possible retraceable paths. The RPA paths result from dressing 
the trivial skeleton-path of length 0. For all other skeleton-paths two types of dressing have 
to be taken into account: 1) additional steps within the skeleton-path (internal dressing), 2) 
irreducible RPA vertex parts inserted at all vertices of the internally dressed skeleton-path. 
To avoid double-counting one has to distinguish between external and internal vertices in 
the internally dressed skeleton diagram. External vertices are connected by one link to 
the skeleton which implies only one forbidden first step for all retraceable-path insertions. 
The dressed "external" vertex C ex (w) is therefore correctly described by Eq. (||). Internal 
vertices, however, are doubly connected to the skeleton-path, such that two first steps are 
forbidden for insertions of the retraceable-path vertex dressings. Therefore, the RPA-dressed 
irreducible internal vertex part is: 



In order to perform the dressing we have to count the number of internal dressings which 
can be considered for a skeleton diagram of length I. For this we define Nf(m,n) as the 
number of paths of m steps on a skeleton diagram of length I, starting at site 0, reaching site 
k (0 < k < I) after step m, and visiting the end-points and I exactly n times (including 
the start at site 0). In terms of these numbers the dressed contribution of a single skeleton 
diagram of length I to the propagator is given by 
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(9) 
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For the numbers iV/^m, n) the following recursive system of equations holds: 



N°(0,n) 
Nf(m + l,n + l) 
Nf(m + l,n) 



Nftm+ l,n+l) 



Ni(m,n), 

N? + \m,n) +Nj t - 1 {m,n), 

(0 < fc < Z) 

iV/ _1 (m,n). 



One can solve these equations by introducing the generating functions 

m,n 

and rewriting the equations (JTTJ) in terms of them. One obtains 



9o{x,y) = V + xygi{x,y), 

9k(x,y) = x(g k -i(x,y) +g k+1 (x,y)), 

(0<k< I), 

gi(x,y) = xygi-x(x,y). 



The general solution of Eq. (|14]) is given by 



9k(x,y) = a+(x,y)r)!l(x) + a-(x,y)r}*(x), 



with 



r]±(x) = ^-(l±VY^Ax^)= V± . 

The coefficients a± are finally determined from the equations for g Q and gi, (|13|) and 
and one obtains 



9i(x,y) 



xy 2 (r}+ - r]-) 



rf + (1 — xyr]_) 2 — rjL (1 — xyi] + ) 
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This generating function can be viewed as the Green's function of a particle propagating 
on a linear chain of length / with hopping amplitudes at the edge bonds different from the 
internal ones [ p~5|j . 

From Eqs. ( |T0D and (|12| ) we see that the contribution of a dressed skeleton diagram of 
length I is determined as 

w \ w C in {w) J 

For the special values of the variables of gi(x,y) appearing in Eq. ([Tj]) the expression 
(|I%|) simplifies considerably since: 1 — xyr] + (x) = 0. We therefore obtain the very simple 
final result: 

Gi(w) = G RPA (w) ( 2t ) , (20) 

\w + ^w 2 -w RPA J 

which holds even for I = 0, although our derivation did not include this case. 

In terms of these dressed skeleton diagram contributions we can now write the propagator 

as: 

oo 

G»(w) = XX-CH, (21) 

1=0 

where K\ ^ is the number of different bare background restoring skeleton-paths of I steps 
between the end-points i and j. On a Neel background only skeleton-paths composed of an 
even number of steps / will contribute, as sites i and j must belong to the same sublattice. 
Trivially, considering that Kf, = 5 it j, the RPA is contained in the propagator of Eq. (|21|). 

Eq. ([H|) represents a resummation of the Nagaoka expansion (|l]) motivated by the merits 
of the RPA. As we will show in the next section, it can be used quite efficiently because 
the number of skeleton-paths of length / is much smaller than the total number of paths of 
the same length. As discussed in the Introduction, the contributions from non-retraceable 
paths get less important the higher the dimension of the lattice. Therefore, we have chosen 
the two non trivial systems of lowest dimension to demonstrate the usefulness of Eq. (pi|): 
the double chain and the square lattice. 



The resummation ( pi]) opens a new chance for a qualified discussion of the band-edge 
problem. Brinkman and Rice || have given an argument which would suggest that the 
spectrum of the correlated hole motion on a Neel background extends to the Nagaoka energy 
wfp = Zt, via exponential band-tails. Sorella et al. |14]] appear to provide evidence in favour 
of this scenario by showing that their spectrum extrapolates to the Nagaoka energy. In our 
opinion, both the Brinkman-Rice argument and the Sorella et al. evidence are not conclusive. 
It is, in fact, obvious from a total spin decomposition of the Neel state that the spectrum 
extends to wfp for finite size systems: the maximum spin S = L/2 is contained in a L-site 
Neel background with a probability of order 2~ L . Therefore the issue is not at all whether the 
spectrum extends to wpp for systems of sufficient size, as demonstrated by the extrapolation 
ITU , but rather whether this part of the spectrum survives with a finite weight in the 
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thermodynamic limit. The results we present in the next section point to band-edges below 
w F p- 

The position of the band-edge is related to the radius of convergence of the series in 
Eq. (pip. If we replace the series by any finite summation the band-edge remains at the 
RPA value (Q). An extension of the spectrum beyond this value can only result from the 
series diverging for energies w < w c with w c > wrpa- One can quite generally assume an 
exponential increase in the number of skeletons with the path-length: 

K\. oc a\ (I -» oo). (22) 

This determines the radius of convergence, which using Eq. (|2^) gives the following simple 
relation between the position of the band-edge w c and the asymptotic growth parameter, a, 
of the number of skeleton-paths: 



, w RPA , a < VZ - 1 
w c = { (23) 



The band-edge energy given by (^3]) grows monotonically with increasing a until it 
reaches wfp for a = Z — 1. The total number of non-retraceable skeleton-paths of length / 
for the hole, irrespective of whether the spin background is restored or not, is Z{Z — 
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Thus a = Z — 1 would mean that the number of Neel background-restoring skeleton-paths 
has the same growth behaviour as the total number of skeleton-paths. This is extremely 
implausible at dimensions D > 3. In this case the majority of non-retraceable closed random 
paths is largely free of points visited more than once, which means that to restore a Neel 
background the hole essentially has to walk twice through such a path. This implies that 
a should be close to \J Z — 1, rather than to Z — 1. The above argument gets better the 
higher the dimensionality. 

We conclude from the above consideration that we should expect the band-edge to be 
close to wrpa-i that i.e. far away from wfp, for lattices at high dimension. At the same 
time, we wish to emphasize that we have not been able to make the above argument exact 
by a strict estimate of an upper bound smaller than Z — 1 for the parameter a. We will see, 
however, in the next section that, even for the low dimensional systems we have considered, 
our numerical data provide good evidence for the conjecture that a < Z — 1 is generally 
true. 

To calculate exactly the propagator given by Eq. (^T|), one would need to determine the 
number of distinct background-restoring bare skeleton-paths of all lengths, which is a very 
complicated problem. In the next section we describe how we have overcome this problem 
for the concrete cases of a double chain and a square lattice Neel backgrounds. 

III. RESULTS AND DISCUSSION 

As mentioned in the previous section, we did find a way to estimate the numbers of 
different bare skeleton-paths for all lengths and thus, through Eq. ([H]), obtain the propagator 
of a hole on a Neel background in two special cases. In this section, we will first describe our 
estimation in general, and in the two following subsections present and discuss the results 
we obtained for the density of states and spectral density, in the cases of a double chain and 
a square lattice. 

To be specific, we first determined numerically the exact numbers of different Neel 
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background-restoring skeleton-paths for as many path-lengths as our computer facilities 
allowed. For the case of a double chain we obtained all skeleton-paths up to a length of 
32 steps, and for lengths 34 and 36 we only determined the closed paths. The results are 
presented in Table |. For the square lattice Neel background we obtained all skeleton-paths 
up to length 24. In Table [IT] we exhibit these numbers. 

From the numbers of skeleton-paths exactly obtained, we determined the asymptotic 
behaviour as a function of path-length exhibited by our data. This we could adequately fit 
by the following functional form depending on three parameters: 

4;-%^ ( 24 ) 

Notice the exponential dependence with the length which leads to the band-edge value given 
by Eq. (0), as described in the previous section. We determined the parameter a from fitting 
the data for the local propagator, the numbers of closed skeleton-paths. The band edge of 
the energy spectrum which we determine from the local propagator is a "global" quantity, in 
the sense that it has to contain all poles of the Green's function, and the non-local Green's 
functions should not extend beyond this edge. Taking the latter fact into account and also 
to avoid spurious singularities in the self-energy obtained, we took the parameter a from the 
local propagator fit and fixed it while fitting the non-local propagators. The denominator in 
Eq. (0), depending on the parameter (3, will essentially account for the energy-dependence 
of the propagator near the band-edge. Finally we included a proportionality constant, the 
parameter C. 

In Tables |[| and |TV] we detail the asymptotes obtained for the double chain and the 
square lattice respectively. The asymptotes were obtained by a weighted fit to Eq. (|24]) of 
our exact short range data, attaching to the latter a dispersion decreasing with path-length. 



We employed a general linear least squares fit solved by singular value decomposition [jTo 
In particular, for the numbers of different skeleton-paths of I steps extending from the origin 
% to the end-point j defined by the propagator considered, which were obtained numerically, 
we considered a standard deviation o\. 
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<j\ oc e 



-0.3/ 



(25) 



Considering all data equally weighted one arrives at asymptotes similar to those presented 
here, but taking a path-length dependent weight attached to the exact short range data will 
certainly produce a better representation of the asymptotic behaviour for long path-lengths. 
To quantify this, we can mention that the relative deviation of the number of skeleton-paths 
from the asymptotes tabulated is of the order of 10 _4 (double chain) and 10 _3 (square lattice) 
for the longest paths exactly investigated. With the weighted fit (p^| ) which we employed we 
could minimize spurious effects in our results, such as negative peaks in the spectral density. 
It is important to remark that the a parameter we obtain, which determines the band-edge, 
is very stable to the consideration of different fits. We checked this by comparing weighted 
fits ( p5[ ) with other exponents of e~ l , including the equal- weight case, and considering other 
functional forms for the standard deviations of the data. In all cases the a parameter 
obtained was well within 2% of the value presented in this work for the double chain, and 
8% for the square lattice where fewer exact data from which to obtain the asymptote are 
available. 

Having obtained the asymptotes for the skeleton-path numbers, and taking into account 
that for short path-lengths (say up to length Z ) we know the exact values of these (K) 
numbers, we evaluated the propagator given by Eq. (^TJ) in the following way: 



The second term in the above expression represents our estimation of the asymptotic or 
long skeleton-path contributions extrapolated to all path-lengths, while the last term is a 
correction to include exactly the short range contributions evaluated numerically. It is an 
important detail that the infinite series appearing in Eq. (|26|) is related to the polylogarithm 
function: 





(26) 



z 



ffc) = E 



(27) 



71=1 
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for which the analytic continuation beyond the radius of convergence (\z\ = 1) is well known. 
We obtained the numerical values for this function employing the Mathematica program [|T7| . 



Having evaluated the local propagator as described by Eq. ( p6|) we determined the density 
of states p(w): 

p(w) = —~Im[Gii(w)], (28) 

7T 

and the Fourier transform of the propagator in real space as: 

GQl,w) = » ,kM: Go.A>c). (29) 
We define the spectral density as: 

A(k,w) = -~Im[G(k,w)], (30) 

7T 

and the self-energy E(k, w) by: 

G(Kw) = (31) 

w — L(k, w) 

In the following we will present the results we obtained for these quantities, and discuss the 
similarities with the results which Sorella et al. [13] obtained employing the Lanczos spectra 
decoding method. 



A. Double Chain 

In Fig. we plot the density of states obtained for the double chain with the method 
described above, and, to compare, include the RPA density of states in the figure. We 
obtain a shift of the band edge, w c = 2.94t, from the RPA value {wrpa = 2.83t) towards the 
Nagaoka energy (wpp = 3t), in accordance with the a parameter obtained (see Table [TTT]): 
Qrpa = v^2 < a ~ 1.88 < app — 2. This result was referred to already in Section [II]. The 
good estimation we obtain for a (independent of the asymptotic fit within 2%) and thereby 
for the band-edge, indicates the absence of band-tails extending to the Nagaoka energy for 
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the double chain. The parameter f3 determines the shape of the density of states near the 
band-edges via p(w) oc (w^ — w 2 ) 13 ^ 1 . 



We obtain a density of states remarkably similar to Sorella et al. [I^j, not only in the 
main features present but quantitatively as well. In particular, there is a broad peak-like 
structure near the band-edge, which we identified as due to the asymptotic behaviour, second 
term in Eq. (|26|) , corresponding to long skeleton-path contributions to G^. In fact, short 
range contributions, which are exactly accounted for by the third term in Eq. (|26|) , are pro- 
ducing a partial cancellation of the peak which would be obtained due to long skeleton-path 
contributions. In the density of states we can also identify a broad shoulder at intermediate 
energies and a second shoulder at lower energies as arising due to the short skeleton-path 
contributions, which we exactly calculate here. 

To compare the orders of magnitude of the different contributions involved in G(k, w) of 
Eq. fl29|), we depict in Fig. [|the imaginary part of the local propagator and the most relevant 
non-local Green's functions between nearest neighbours. We recognize a quick decrease in 
relative weight as one considers neighbours further apart. This means that we obtain good 
convergence for the Fourier transform of the Green's function, G(k, w), by including the 
contributions up to tenth nearest-neighbours. 

In Fig. H] we plot the spectral density ( pOf) as a function of energy, for three different 
values of momentum k. The spectral weight distribution is in general agreement with the 



results of Sorella et al. |14|]. Concretely, for k = (0,0) we obtain a peak-like structure near 



the band edge and a broad shoulder at lower energies, both features present in Ref. |14| with 
slight differences in magnitude. We investigated the origin of those features, and found that 
the band-edge peak is determined by the asymptotic or longer skeleton-path contributions 
to the k-dependent Green's function (where the local and the non-local propagator to 2 nd 
nearest neighbours are the dominant terms). The short range contributions tend to narrow 
and decrease the weight of the peak. It's interesting to note that the non-local propagators 
contributing to G(k, w) tend to shift the peak arising from the local propagator towards 
the band-edges and to increase its weight. Meanwhile, the broad shoulder at intermediate 
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energies originates from the short range contributions to the non-local propagators (mainly 
those corresponding to 2 nd and 3 rd - nearest neighbours). In our case the main peak appears 
accompanied by a narrow secondary peak which is a totally spurious trace of the RPA band 
edge, mainly appearing in the non-local propagators, of no relevance whatsoever. 

To exhibit the dispersion, in Fig. [3] we also plot the spectral density for k = (^-,0) and 
k = (7r,0). Again here we observe the same weight distribution exhibited by the data of 
Sorella et al. [f[4j| , where A(k = (^-,0),w) contains the highest peak-like structure near the 
band-edge in the whole Brillouin zone, while this structure is broadest for the k = (tt, 0) 
data. Also the peak-shifts towards the center of the band agree with those obtained in 



Ref . JTJ] . For the k ^ cases, we find that the short range part of the non-local propagators 
contributing to G(k, w) are increasing the weight of the peak which would arise from the 
local propagator and shifting it towards the center of the band. The same short range part 
of the non-local propagators is responsible for the decrease in spectral weight obtained at 
intermediate energies (where for k = (0, 0) the broad shoulder appeared) and for the increase 
of weight in the center of the band. 

Finally, in Fig. [| we plot the self-energy obtained for the same k-values taken for the 
spectral density results, and, for purposes of comparison, we include the dispersionless RPA 
self-energy. The effects of dispersion are clearly observable, and important departures from 
RPA are obtained. In Fig. |] (a) we plot the real part of the self-energy. The k = (0, 0) 
data show maximum departures from RPA for energies near the center of the band, while 
for the other k-values these appear shifted towards intermediate energies. The origin of 
these departures lies mainly in the short range parts of the non-local propagators. Similar 
comments apply to the imaginary part of the self-energy plotted in Fig. |] (b). For k = (0, 0), 
we obtain an important reduction of the lifetime as compared to the RPA for energies near 
the center of the band. 



16 



B. Square Lattice 



In Fig. |5] we plot the density of states we obtained for the square lattice. As for the double 
chain, we obtain a shift of the band-edge, w c = 3.61t, from the RPA value {wrpa = 3.46t) 
towards the Nagaoka energy (wfp = 4t). This is determined by the a parameter obtained 
from the asymptotic fit (see Table [TV]) through Eq. (^): a RPA = \f?> < a ~ 2.31 < app = 4. 
As commented before, the reliability of our estimation for a (within 8%) clearly indicates 
the absence of band-tails extending to the Nagaoka energy in the spectrum of the square 
lattice. This result also confirms our conjecture of Section [TI| that the band-edge departures 
from the Nagaoka energy increase with dimensionality. 

The almost featureless density of states, with minor departures from the RPA result as 
expected from the estimations in terms of the dimension mentioned in the Introduction, 
agrees with that obtained by Sorella et al. |14| apart from some minor oscillations. These 
oscillations result from the short skeleton-path contributions to the local propagator. In 
contrast to the double chain, here the only observable departure from the RPA density of 
states resulting from the asymptotic behaviour is the band-edge shift. 

In Fig. ^ we depict the imaginary part of the local propagator, and the most relevant 
non-local Green's functions between nearest-neighbours. Again, a quick decrease in weight 
appears on considering neighbours further apart, but this plot shows that non-local prop- 
agators have much less weight than the local one in comparison to the double chain (see 
also Fig. |H). This agrees with the estimation of the importance of non-local corrections to 
the RPA propagator as a function of the dimension mentioned in our Introduction. For the 
square lattice we obtained the spectral density including in the Fourier transform of the 
Green's function, G(k, w), contributions up to those of twelfth nearest-neighbours. 

In Fig. [7] we plot the spectral density as a function of energy, for three different values of 
momentum k. Here for k = (0, 0) we obtain a low peak-like structure near the band edge and 
a shoulder at lower energies, which is broader than the one in Ref. |14| where the data seem 
to have been insufficient to enable them to trace the complete spectral density curve. Both 
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structures result mainly from the short range part of the non-local propagators contributing 
to G(k,w). The effect of dispersion on the spectral weight distribution is coincident with 
that shown by Sorella et al. |TJj]. In Fig. |7| we also plot the spectral density for k = (n, 0), 
which exhibits the highest and broadest peak-like structure, and k = (|, ~), which shows 
three maxima, all in accordance with Ref . |14[ . Again, the short range part of the non-local 
propagators contributing to G(k, w) accounts for this distribution of weight. 

In Fig. H we plot the self-energy obtained for the same k-values taken for the spectral 
density, including the dispersionless RPA self-energy for hole motion on a square lattice with 
Neel spin order. The qualitative effects of dispersion are similar to those obtained for the 
double chain, but the departures from RPA obtained are smaller, as expected. As before, the 
departures obtained mainly result from the short range part of the non-local propagators. 

IV. SUMMARY 

In this work we studied the motion of a hole on a Neel background, neglecting spin 
fluctuations, in the framework of the Nagaoka expansion for the Green's function. Starting 
from the retraceable-path approximation, known to become exact in the one- and the infinite 
dimensional cases, where no dispersion appears, we derived a resummation of the Nagaoka 
expansion by considering non-retraceable skeleton-paths and dressing them by retraceable- 
path insertions. The contribution of each of such dressed skeleton-paths was evaluated 
exactly. The problem is then reduced to the determination of the numbers of different bare 
skeleton-paths of each length. We numerically exactly obtained these numbers up to length 
36 for the double chain, and 24 for the square lattice. From these numbers, we determined 
an asymptotic extrapolation for all path-lengths, and used all this information to evaluate 
the Green's function. We then determined the density of states and spectral density for 
the double chain and square lattice cases. We determine the band-edges through a growth 
parameter, a, obtained in the asymptotic fit. The reliability of our determination of this 
parameter, placing the band-edge at a value between the RPA edge and the Nagaoka energy, 
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points towards the absence of band-tails extending to the Nagaoka energy in the spectrums of 
the double chain and the square lattice. At the same time, this confirms our conjecture about 
the increase of the band-edge departures from the Nagaoka energy with dimensionality. Our 
results deviate from the exact solution of the problem only due to the differences between 
our asymptotic extrapolation and the unknown exact numbers of skeleton-paths longer than 
36 steps for the double chain and 24 for the square lattice. That these deviations are small 
is confirmed by the general coincidence of our density of states and spectral density results 
with those obtained through the Lanczos approach by Sorella et al. [|14]]. Our "analytic" 
approach has the added advantage of determining the relevant hole-paths responsible for 
the main features present in the density of states and the spectral density. 
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FIGURES 

FIG. 1. Double chain. Density of states: our result (full line); retraceable-path approximation, 
RPA, (dot-dashed line). 

FIG. 2. Double chain. — Imaginary part of the local propagator (full line), of the RPA(local) 
propagator (dot-dashed line); of the non-local propagator between the origin and 2 nd (dashed line), 
3 rd (boxes), and 6 th nearest neighbours (crosses). 

FIG. 3. Double chain. Spectral density A(k, w) as a function of energy for: k = (0, 0) (full 
line), k = (^,0) (dashed line), k = (tt,0) (dot-dashed line). 

FIG. 4. Double chain. Self-energy X(k, w) of the hole as a function of energy: (a) Real part, 
(b) - Imaginary part. RPA (boxes); k = (0,0) (full line); k = (2f,0) (dashed line); k = (vr,0) 
(dot-dashed line). 

FIG. 5. Square lattice. Density of states: our result (full line); retraceable-path approximation, 
RPA, (dot-dashed line). 

FIG. 6. Square lattice. — Imaginary part of the local propagator (full line), of the RPA(local) 
propagator (dot-dashed line); of the non-local propagator between the origin and 2 nd (dashed line), 
3 rd (boxes), and 5 th nearest neighbours (crosses). 

FIG. 7. Square lattice. Spectral density A(k, w) as a function of energy for: k = (0, 0) (full 
line), k = (n, 0) (dashed line), k = (|, |) (dot-dashed line). 

FIG. 8. Square lattice. Self-energy S(k, w) of the hole as a function of energy: (a) Real part, 
(b) - Imaginary part. RPA (boxes); k = (0,0) (full line); k = (ir,0) (dashed line); k = (f , f ) 
(dot-dashed line). 
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TABLES 

TABLE I. Exact numbers of non-retraceable Neel background-restoring skeleton-paths on a 
double chain, classified according to path- length or number of steps (1), and distance (d) between 
origin and end-point, measured in units of the lattice parameter. Due to symmetry in the double 
chain, n s d2 = 2 different sites are end-points with the same 0), thereby contributing equally to 
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TABLE II. Exact numbers of non-retraceable Neel background-restoring skeleton-paths on a 
square lattice, classified according to path- length or number of steps (1), and distance (d) between 
origin and end-point, measured in units of the lattice parameter. Due to symmetry, n s d2 different 
sites are end-points with the same d (see row in brackets). 
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TABLE III. Double chain: Asymptotes obtained for K\ j (as denned by Eq. (|24|)). In the first 
column we enter the relative position of end-point j with respect to the origin i, denoting by n th (i) 
a n th nearest-neighbour site of the origin i. For clarity, in the second column we write the distance 
d between i and j (squared). The last column indicates which (n l d2 /n s d2 ) data (from Table |) were 
employed to determine each asymptote by the weighted fit (Eq. (p5|)), by stating the respective I. 



j d 2 a (3ij Qj I 

j = i 1.88136 2.39084 1.26191 22-36 

2 nd (i) 2 " 2.27482 0.66616 22-32 

3 rd {i) 4 " 2.17748 0.33958 22-32 

6 th {i) 10 " 1.76242 5.278xl0" 2 22-32 

7 th {i) 16 " 1.42770 9.05xl0~ 3 22-32 

10 th (i) 26 " 0.86016 5.7xl0" 4 22-32 
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TABLE IV. Square lattice: Asymptotes obtained for Kjj (as defined by Eq. (|24|)). In the first 
column we enter the relative position of end-point j with respect to the origin i, denoting by n th (i) 
a n th nearest- neighbour site of the origin i. For clarity, in the second column we write the distance 
d between i and j (squared). The last column indicates which (n d2 /n d2 ) data (from Table ||) were 
employed to determine each asymptote by the weighted fit (Eq. (p5|)), by stating the respective I. 
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